Positron-molecule annihilation by capture into vibrational Feshbach resonances of 

infrared-active modes 
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Enhanced positron annihilation on polyatomic molecules is a long-standing and complex problem. 
We report the results of calculations of resonant positron annihilation on methyl halides. A free 
parameter of our theory is the positron binding energy. A comparison with energy-resolved anni- 
hilation rates measured for CH 3 F, CH 3 C1, CH 3 Br [Barnes et al. Phys. Rev. A 74, 012706 (2006)] 
shows good agreement and yields estimates of the binding energies. 
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In this paper we calculate the positron-molecule an- 
nihilation rate due to resonant capture of positrons 
by infrared-active vibrational modes, and observe good 
agreement with recent experimental data for methyl 
halides |l(. 

When a fast positron interacts with matter, it under- 
goes a quick succession of ionizing and other inelastic col- 
lisions, and slows down to eV or thermal energies before 
annihilation. The low-energy positron annihilation rate 
in a gas with number density n is usually parametrized 



A : 



irr cnZ cfi , 



(1) 



where a a is the annihilation cross section, v is the 
positron velocity, c is the speed of light, r is the clas- 
sical electron radius, and Z c s is an effective number of 
electrons per gas atom or molecule, that contribute to 
annihilation [2j, |3j ■ Originally, Z c g was introduced in ex- 
pectation that the annihilation rate would be in propor- 
tion to the number of target electrons, Z. However, early 
experiments 0, 0, @ and later systematic studies 0, IS El 
found that for many polyatomic molecules Z cS exceeded 
Z by orders of magnitude. It also showed strong chemi- 
cal sensitivity and rapid growth with molecular size (see 
review [1(3]'). 

Explanations of high molecular Z e R were sought 
in terms of positron virtual or weakly bound states 
|ll| , resonances 0, , long-lived vibrationally excited 
positron- molecule complexes 0, and virtual Ps forma- 
tion [l^ . At the same time, annihilation calculations 
which neglected molecular vibrations, failed to reproduce 
"anomalous" Z c $ for polyatomics [la, [lj| , but gave evi- 
dence that Z c ff depend on the molecular geometry |l7j . 

These efforts highlight the fact that positron-molecule 
annihilation is a complex problem. Nevertheless, a the- 
ory developed in Refs. 0, 0, provides a frame- 
work for analyzing this phenomenon. There are two ba- 
sic mechanisms of positron annihilation, direct and res- 
onant. The direct mechanism applies to both atoms 
and molecules and involves annihilation of an incident 
positron "in flight" . Its contribution is enhanced when 



a low-lying virtual or weakly bound positron state is 
present, leading to Z e g up to 10 3 for room-temperature 
positrons [HIi^. 

Resonant annihilation occurs for molecules capable of 
binding the positron. To be captured into a bound state, 
the positron energy must be absorbed by a vibrational 
excitation of the positron-molecule complex. This gives 
rise to a vibrational Feshbach resonance (VFR) at the 
incident positron energy e = E v + Eq, where E v is the 
vibrational excitation energy, and £o < is the positron 
bound state energy. The positron bound in the VFR can 
annihilate (or undergo detachment). The probability of 
annihilation is proportional to the resonance lifetime. 

For non-monoenergetic positrons and closely spaced 
resonances, their contribution to Z e g is proportional to 
the vibrational level density at E ss e — £o 0, 0] . If 
the positron VFR were due to excitation of fundamentals 
alone, this density would be proportional to the number 
of modes. Experimental Z Q g show much faster increase 
(e.g., Z cfi = 3500, 11300, and 37 800 for C 3 H 8 , C 4 Hi , 
and C5H12, respectively). This means that positron at- 
tachment involves excitation of overtones and combina- 
tion vibrations. Large Z e f{ are then related to high total 
vibrational spectrum densities in the polyatomics. 

The important role of vibrations was recently verified 
by measuring the energy dependence of Z e g at sub-eV 
energies with a high-resolution positron beam 0, Eij ], 
These experiments uncovered peaks in Z c ff, whose en- 
ergies corresponded to those of molecular vibrational 
modes. In particular, for all alkanes larger than methane, 
Z e g displayed a prominent C-H maximum. Its down- 
shift from the C-H mode energy (0.37 eV) provided a 
measure of the positron binding energy 23]. Observa- 
tion of such peaks means that excited fundamentals act 
as vibrational doorway states |20| , leading to multimode 
vibrations through intramolecular vibrational relaxation 
(IVR). 

Therefore, to compute Z c ff for polyatomics, one must 
account for strong electron-positron correlations and 
positron binding, the interaction between positronic and 
vibrational degrees of freedom and intramolecular vibra- 
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tional mixing. This makes ab initio calculations of high 
molecular Z e g very difficult. However, as we show below, 
for small polyatomics some basic features of resonant an- 
nihilation can be tested by relatively simple calculations. 

The resonant part of the annihilation cross section can 
be written using the Breit-Wigner formalism 

13 HSU, 



-Y 



a V a Y e 



(e - E v 



so) 



•rj/4' 



(2) 



where r", T e w , and T u are the annihilation, elastic and 
total widths of vth resonance, g v is its degeneracy, and k 
is the positron momentum (atomic units are used). The 
annihilation width is proportional to the electron density 
at the positron in the positron bound state, p ep , 

T° = irrlcp ep . (3) 

From Eqs. and (J2J, the resonant Z e s is given by 



7 (rcs) _ 7T >p 



9^1 



{e-E u -e ) 2 +Tl/4 



(4) 



We will now use this equation to calculate the contribu- 
tion of infrared-active modes to B ' . 

Consider a compact polyatomic molecule that can bind 
the positron with a small binding energy \eq\ = k 2 /2 <C 1 
eV. The wavefunction of the bound positron is very dif- 
fuse and behaves as ipo — Ar~ 1 e~ Kr outside the molecule. 
Since large distances dominate, the normalization con- 
stant is given by A ~ {k/2tt) 1 I 2 |2^ . 

Suppose that the vibrational modes in this small-sized 
polyatomic are not mixed with overtones or combination 
vibrations. Given the smallness of the binding energy, the 
vibrational excitation energies of the positron-molecule 
complex should be close to the fundamental frequencies 
lo v of the neutral molecule. In this case the sum in Eq. 
(@J) is over the modes is, and E v ss lo v . Some (or even 
all) of these modes can be infrared active. The positron 
capture into such excited states is mediated by the long- 
range dipole coupling. This allows one to calculate their 
contribution to Z^ s ^ . 

Consider a positron with momentum k incident on the 
molecule in the vibrational ground state <I>o(R), where 
R represents all the molecular coordinates. If k 2 /2 as 
u>i/-\-£q, the positron can be captured in the VFR, where 
it is bound to the molecule in a vibrationally excited state 
$„(R). The corresponding width Tf, can be found from 



ri = 2n / \A vk \ J 8{k*/2-w v - e ) 



d 3 k 
(2tt) 3 ' 



(5) 



where A v ^ is the capture amplitude. We calculate it by 
using a method similar to the Born-dipole approximation 



|2a|2J, as 



4iri d,, • k 



dr 



r $o(R)drdR 



1 5 k 2 
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(6) 



where d is the dipole moment operator for the molecule, 
= ($ w |d|$o), and 2F1 is the hypergeometric function 
Substitution of Eq. @J into Eq. © gives 



(J) 



where h{£) = C 3/2 (l - 0" 1/2 [2^ (|, 1; §5 - 0)] * 

is a dimensionless function of £ = 1 + Eq/uj^, such that 
£(0) = £(1) = 0, and Vax ~ 0.75 at £ w 0.89. 

Equation Q shows that the elastic width of a positron 
VFR for an infrared active mode is basically determined 
by its frequency 10 v and transition dipole amplitude d u , 
known from infrared absorption measurements [29j. 

For weakly bound positron states the density p ep is a 
linear function of k [ljj . It can be estimated as 



Pep = (F/2tt)k, 



(8) 



with F « 0.66 0| . The same constant characterizes the 



contribution of direct annihilation, Z. 



(dir) 



|l9l | . It is enhanced at small positron momenta by 



F/(k 2 + k 2 ) 
;nta by the 

presence of a weakly-bound (or virtual) state [lj, yfj . 

In a recent paper Q measurements of Z c ff for CH3CI 
and CH^Br using a cold trap-based positron beam, have 
been reported. The energy dependence of Z e g for these 
molecules (and CH3F measured earlier j^) shows peaks 
close to the vibrational mode energies. This points to 
an important contribution of resonant annihilation in all 
three molecules, although the maximum Z e g value for 
CH3F (250) is much lower than those for CH3CI and 
CH 3 Br (1600 and 2000, respectively). 

These molecules have Cs v symmetry, and all six of 
their vibrational modes are infrared active (see Table|l]for 
CH3CI). Methyl halides are also relati vely small, which 
means that IVR may not take place |31|. This makes 
them ideal for application of our theory. Equations 
ijjj , and (JHl allow one to calculate the contribution of all 
VFR to Z^fi , and the only free parameter of the the- 
ory, i.e., the positron binding energy, can be chosen by 
comparison with experimental Z e g. 

In order to do this, Z^ s ' from Eq. Q must be aver- 
aged over the energy distribution of the positron beam 
|22| . The latter can be modelled by a combination of the 
Gaussian distribution in the longitudinal direction (z) 
and Maxwellian distribution in the transversal direction 
(_L). The corresponding probability density of the total 
positron energy, s± + s z , is 



fe{£±,£z 
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kBT±V 2na 2 



exp 



£_l 
k B T± 
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TABLE I: Characteristics of the vibrational modes of CH3CI. 



Mode 


Symmetry 


9u 


a 


d v 










(meV) 


(a.u.) 


(a.u.) 


V\ 


ai 


1 


363 


0.0191 


4.87 x 10" 6 




ai 


1 


168 


0.0176 


1.91 x 10~ 6 


vz 


ai 


1 


91 


0.0442 


6.52 x 10~ 6 


"i 


e 


2 


373 


0.0099 


1.34 x 10" 6 




c 


2 


180 


0.0162 


1.74 x 10" 6 


V6 


c 


2 


126 


0.0111 


5.66 x 10~ 7 



"Mode energies u„ and dipole amplitudes d„ from Ref. I'Jl.l . 

where k B is the Boltzmann constant, T±_ is an effective 
transversal temperature of the beam, e is the mean lon- 
gitudinal energy of the positrons, as measured by the 
retarding potential analyzer, and a = (5 z /- v / 81n2, 8 Z be- 
ing the full width at half-maximum. The values of k B T± 
and S z are taken from experiment to be 25 meV. 

The averaging, Z^ s) (e) = J Z^ s) f e (sj_,e z )d,£±de z , 
can be done analytically, since the widths of the reso- 
nances, T^ = Y e v + T", are small compared to the en- 
ergy spread of the positron beam. Indeed, values from 
the last column of Table [i] show that the elastic widths 
of the VFR are less then 0.1 meV. Typical annihilation 
widths are even smaller. For example, for a binding en- 
ergy of 10 meV (k = 0.027 a.u.), Eqs. © and © yield 
r° = 3 x 10~ 9 a.u. = 0.1 jtxeV. These estimates also show 
that r„ w T e y: i.e., that the total width of the resonance 
is dominated by its elastic width. 

Hence, to integrate over e± and e z we replace the Breit- 
Wigner profiles in Eq. by <5-functions, and obtain 



^eff 



(e) 



(9) 



where £„ = fc^/2 = uj v + e$ is the resonance energy, and 



A(E) = 



1 



k B T_ 
1 + $ 



exp 



2{k B T 1 _) 
— (-I 



exp 



k B Tj_ 



E 



k B T± 



(10) 



with $(2;) being the standard error function. 

The function A(E) is a convolution of the 5-function 
with the positron energy distribution. It describes the 
appearance of a narrow resonance when measured with 
the trap-based positron beam, and is shown in Fig. ^ 
Due to the transversal energy component, its maximum is 
downshifted by 12 meV from the true resonance position. 
The shape of A(E) is also markedly asymmetric, with 
an extended low energy tail. It agrees well with those of 
the observed C-H peaks 0)H3|. Note that the positron 
energy distribution was taken into account in experiment 
by assuming a 16 meV difference between the positron 
total and longitudinal energies |22^ . 
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FIG. 1: Comparison of the resonance shape function A(E), 
for fc^Tx = 8 Z = 25 me V (curve) with the measured C-H 
peak in propane (circles) pH |22^| . For comparison, experi- 
mental Z a g has been scaled vertically and shifted horizontally. 
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In Fig. |21 we compare the beam-energy-averaged the- 
oretical Z e s from Eq. @ added to the direct contribu- 
tion Z^ T ^ , with measured Z e g for methyl halides Q, • 
Theoretical curves have been obtained using the binding 
energy of \e a \ = 0.3, 25, and 40 meV, for CH 3 F, CH 3 C1 
and CH 3 Br, respectively. 
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FIG. 2: Comparison between experimental Z e g (•, Ref. Q]) 
and theoretical Z e a obtained as a sum of the direct and res- 
onant contributions, using the binding energy Eo — 0.3 meV 
(CH 3 F), 25 meV (CH3CI), and 40 meV (CH 3 Br). Vertical 
bars show the energies of molecular fundamentals. 
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Given the complexity of the problem and the fact that 
eg is the only free parameter in the calculation, the agree- 
ment between theory and experiment in Fig. [21 is re- 
markable. In accord with Eq. @, every vibrational 
mode gives rise to a VFR, whose relative magnitude is 
determined by the factor g u /k u (since T^/r^ « 1). On 
the positron longitudinal energy scale, the resonances are 
downshifted from the mode energies by the positron bind- 
ing energy and a further 12 meV due to the positron 
energy distribution. 

Besides determining the resonance positions, the bind- 
ing energy also affects the overall magnitude of Z^£ s ' via 



Pep OC |£q 



1/2 



[Eq. 



Hence, the smallness of Z c ff in 



CH3F in comparison with those of CH3CI and CH^Br is 
related to the weakness of its binding. This is in turn 
related to the smaller dipole polarizability and higher 
ionization potential of fluoroniethane, which make it less 
attractive for the positron. 

Note that the infrared absorption strengths of the 
modes and the corresponding elastic widths, T e y ~ 0J v d1, 
may vary considerably from mode to mode. On the other 
hand, the contribution of different modes to Z c r are simi- 
lar, apart from energy shift and g v lk v factor. As a result, 
the energy dependence of Z e g has little resemblance to 
the molecular infrared absorption spectra [l|. The rela- 
tion r„ w Y e v also means that the contributions of the 
VFR are not sensitive to the exact values of the elastic 
widths. Therefore, our use of the "Born-dipole" approx- 
imation in the derivation of Eq. JJJ is not expected to 
lead to sizeable errors in ZjS*^ . 

In conclusion, we have presented a theory of positron 
annihilation by capture into vibrational resonances of 
infrared-active modes. It agrees well with measured Z Q g 
for methyl halides and yields estimates of the positron 
binding energies for these molecules. 

This theory can also be used to investigate the con- 
tribution of infrared-active-mode VFRs to Z e g in other 
small polyatomics that can bind positrons. Such calcula- 
tions will likely underestimate the Z e s, because the res- 
onances associated with other (non-dipole) modes may 
contribute just as much, as long as their elastic widths 
are greater than the annihilation width. 

In molecules where multiquantum vibrations are cou- 
pled by anharmonicity, the number of VFRs populated 
by positron capture will be greatly increased, leading to 
much higher Z c g. However, the same coupling will also 
allow the VFR to decay by positron emission to vibra- 
tionally excited states of the molecule. This will increase 
the total resonance widths, thereby reducing their indi- 
vidual contributions. Calculation of Z e f{ for molecules 
with IVR is the next big challenge for the theory. 

The authors are grateful to C. M. Surko and 
J. A. Young for a most helpful discussion of the positron 
energy distribution and valuable comments, and to 
A. V. Korol for discussions. 
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